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Abstract 

A systematic theoretical study of stoichiometric clusters (Nal) n 
up to n = 15 is performed using the ab initio Perturbed-Ion (PI) 
model. The structures obtained are compared to previous pair po- 
tential results, and observed differences between (Nal) n clusters and 
previous ab initio results for other alkali halide clusters are discussed. 
(Nal) n clusters with n up to 15 do not show yet a marked preference 
for geometries which are fragments of the bulk lattice. Instead, stacks 
of hexagonal rings or more open structures are obtained as ground 
structures in clusters with n = 3, 6, 7, 9, 10, 12, 13 and 15, indicat- 
ing that convergence to bulk structure is not achieved yet at this size 
range. Low lying isomers which are fragments of the crystal lattice 
exist, nevertheless, for those cases. The binding energies show that 
clusters with n = (4), 6, 9 and 12 molecules are specially stable. The 
binding energy has been decomposed in contributions which allow for 
an intuitive interpretation. Some electronic properties like ionization 
potentials and electronic energy levels are also studied. 

PACS: 36.40. +d; 61.50.Lt; 71.20.Fi 
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1 Introduction 



The structural and electronic properties of a material depend primarily on 
the state of aggregation]]]]. Small clusters often present significant physical 
and chemical differences with respect to the molecule and the bulk material. 
Cluster studies can provide valuable insight on the development of the prop- 
erties of condensed matter from those of isolated atoms or molecules. One 
important source of information about cluster properties is provided by the 
abundance patterns emerging from the mass spectra. The relative intensities 
of cluster peaks are interpreted as being indicative of variations in the cluster 
stability with size. The measured population usually shows large drops after 
some particular sizes (often the population is also a maximum at those sizes) 
and those particularly stable clusters are known as "magic clusters". The 
magic numbers have been explained by using either electronic or structural 
models. For example, Van der Waals clusters X n appear to form compact 
icosahedral structures which exhibit enhanced stabilities for completed lay- 
ers at values of n=13,55,147,... . In contrast, the stability variations of alkali 
metal clusters have been interpreted in terms of the electronic-shell model|3|]. 
A direct experimental determination of the spatial distribution of atoms or 
molecules in a cluster is a very difficult task; however the knowledge of the 
cluster structure is a prerequisite for a quantitative discussion of many prop- 
erties. Drift tube studies, which measure the mobility of a cluster ion through 
an inert buffer gas under the influence of a weak electric field, provide valu- 
able information about the cluster geometries^], 0. A theoretical study of 
the most stable structures and their evolution as a function of the number of 
atoms or molecules in the cluster is thus of the upmost importance for the 
interpretation of the mobility experiments. 
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There are several experimental techniques for the production of small clus- 
ters of alkali- halide materials: particle sputtering||, inert-gas condensation^] 
and laser vaporization [§, H|. I n this paper we are interested in stoichiomet- 
ric (Nal) n clusters. Large non-stoichiometric "ionized" (NaI) n Na + clus- 
ters have been produced and their mass spectrum analyzed by Martin and 
coworkers ||10||. The spectrum is characterized by a periodic structure which 
has been interpreted as revealing a definite pattern of growth: those clusters 
grow as cuboids (clusters with rectangular faces and rt\ x ri2 x 773 atoms, 
where n\, 712 and 773 are the number of atoms along the three edges), and the 
observed intervals between minima in the spectrum correspond to the num- 
ber atoms necessary to cover one face of the cuboid. The mass spectrometry 
techniques detect only ionized clusters, although one can expect that the 
formation pattern of neutral (Nal) n clusters will be analogous (neutral sto- 
ichiometric TT-i x ri2 x 77.3 cuboids are, however, not possible if ni, ri2 and 



773 are all odd numbers). Kreisle and coworkers ||11|| have performed mobility 
measurements for (NaI) n Na + clusters up to 77 = 18. 

Dieffenbach and Martin have calculated the lowest energy structures of 



Nal clusters using pair potential models |12[. In pair potential studies for 
NaCl clusters, Phillips et al. |13| found very small differences in binding 
energy between isomeric configurations of the same cluster. In order to check 
if the approximations inherent in a pair potential formulation have influence 
on the predicted geometries we use the ab initio Perturbed Ion (PI) model, 
in which the radial relaxation of the wave functions of the Na + and I~ ions 
due to the interaction with the environment is taken into account. Here 
we study (Nal) n clusters up to n = 15. To our knowledge, these are the 
first ab initio calculations on such clusters. The PI model, which has been 



applied to study (NaCl) n clusters in previous works|l4|, [15], [16]], had been 
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originally developed for the study of ionic materials in the crystal phase [^7 



The model is based on the Theory of Electronic Separability (TES)|Tq, |19 



and the restricted Hartree-Fock approximation. The quantum-mechanical 
interactions between the ions are treated in the following way: since the 
model assumes from the start that the cluster is built from closed shell ions, 
when those ions approach and the overlap becomes important, the radial 
deformation of the charge clouds is taken into account in a self-consistent 
way. The directional covalent effects are considered less important than the 
radial deformations of the closed shell ions, and this restricts the class of 
systems that can be treated. 

It is also interesting to compare our results with the results for NaCl 
clusters previously studied by the PI model. We want to see the influence 
of the ion size: in the present case the size difference between the I~ and 
Na + ions is much larger than for NaCl clusters. The paper is organized as 
follows. The PI model is presented in section II. The results for the structures 
and binding energies of (Nal) n clusters up to n= 15 are discussed in section 
III, and the conclusions are summarized in section IV. 



2 Perturbed Ion model for clusters 

According to the Theory of Electronic Separability (TES)[|H|, when the 
system under study is composed of weakly interacting groups, its wave func- 
tion can be expressed as an antisymmetrized product of group wave functions. 
If these satisfy strong-orthogonality conditions^, |2T|, the total energy is the 
sum of intragroup, or net, energies and intergroup interaction energies. For 
a cluster with a fixed three-dimensional arrangement of the ions, we consider 
each ion as a different group. All contributions to the total cluster energy 
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from a given group, the active ion A, can be collected in the effective energy 
of this group 

E eff = E net + E E int = E net + E inf i 1 ) 

Here E^ et is the intragroup energy and E^ t is the energy due to the inter- 
action of the ion A with all the other ions. Group wave functions can be 
obtained by minimizing their effective energies if strong orthogonality condi- 
tions are satisfied among the groups. The sum of the effective energies is not 
equal to the total energy of the system because this sum counts twice the 
intergroup energies. However, we can define an additive energy of the ion A 
by: 

A A 1 AS A 1 A 

and then the total energy of the system can be written 

N 

E = E E tm- (3) 

A=l 

The effective energies can be expressed as expectation values of appro- 
priate effective Hamiltonians. Each ion of interest (the active ion A) will be 
a group in the sense of TES. Its wave function can be described by a single 
Slater determinant and this can be obtained by minimizing the expectation 
value 

E? ff (elec) =< i>A I Htff I i>A > (4) 



of the effective electronic hamiltonian 18 



N A 

Hf ff = Y,tiff®+ E r«\ (5) 

i=l ^<3<i<N A 

hif f (i) = T{i) - Z A r-} + £ + P 8 ®], (6) 

9&A) 
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where i, j run over the Na electrons of ion A and S counts the rest of the 
ions of the cluster. T is the kinetic energy operator and Z A is the nuclear 
charge of ion A. Vfjj(i) represents the potential energy of the ith electron (of 
ion A) in the field created by ion S 

K S ff (i) = -Z s rT s } + V g{%) + V$(t). (7) 

The different parts of Vf^ are the electron-nucleus, classical electron-electron 
and exchange parts of the potential energy. The strong orthogonality between 
the orbitals of the active ion A and those of the other ions S is included in 
H e // by means of the projection operator P s (i). For closed shell ions, this 
operator can be written in the form [P^| 

J*(i)=£l#>(-2e?)<#l, (8) 

g€S 

where g runs over all occupied orbitals 4> g of the ion S with orbital energies 
e s . 

9 



We now consider the explicit form of the effective potential V^[17|. For 
the closed-shell ions considered here the classical electrostatic part of this 
potential is given by 

Vg (n) = / p S (r 2 )r^dr 2 (9) 
where p s {r) is the electron density of ion S. The exchange operator can be 



written as the nondiagonal spectral resolution |p2 

Vjfri) =-£ £ El alrn,S> A(l,ab,S) < blm,S\, (10) 

l m=—l a,b 

where | aim, S > are products of spherical harmonics Y™ and primitive 
radial functions for the S ion, and a and b run over the Slater-type orbitals 
of /-symmetry. A(l, ab, S) are the elements of the matrix 

A = S~ 1 ICS- 1 (11) 
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and S and K, are the overlap and the exchange matrices for the S ion in the 
{| aim, S >} basis. More details can be found in the original paper by Luana 
and Pueyo [fL7|l . 

Finally the effective energy of the active ion A is the sum of the electronic 
energy of eq. (4) plus the nuclear term 

E^ ff (nucl) = - £ Z A v s {\ R A - R s |) (12) 

S&A) 

where v (| Ra — Rs I) is the electrostatic potential of the S group at the 
nucleus of the active A ion. This potential contains nuclear and electronic 
contributions. 

The practical method that we have used to apply this formulation to 
alkali-halide clusters is the following one: for a given distribution of the ions 
forming the cluster, we consider one of them as the active ion A (for instance, 
a particular anion and solve the Self- Consistent-Field equations for anion 
A in the field of the remaining ions, which are considered frozen at this stage. 
The solution obtained is transferred to all the anions equivalent to anion A, 
that is, to the anions which have equivalent positions in the cluster. We 
take then a non equivalent I~ anion (anion B) as the active ion and repeat 
the same process. Evidently, since anions B are not equivalent to anions A, 
the energy eigenvalues and wave functions of electrons in anions B can be 
different from those of anions A. We continue this process in the same way 
until all the inequivalent anions have been exhausted. The same procedure is 
then followed for the cations Na + . The process just described is a PI cycle. 
We iterate the PI cycles until convergence in the total energy of the cluster 
is achieved. Although for bulk Nal there are only two types of active ions 
(Na + and this number is larger for (Nal) n clusters. The reason is that, 
in the cluster, Na + (and I~) ions are in several different environments. Basis 
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sets from Clementi and Roetti[[^| for the ions (5s4p(iVa + ) and lls9p5d(/~)) 
have been used in the calculations. We have checked that this election is 
the best, between all the basis sets reported by Clementi and Roetti, by 
carrying out some exploratory calculations for the molecule and for (NaI) G 
with different basis sets. More precisely, we have tested the use of the basis 
sets optimized for Na, Na + , I and I~ . Between the four different possible 
combinations, the basis set for Na + ,I~ led to the lowest total energies for 
the two clusters studied. Inclusion of diffuse functions did not show any 
substantial improvement. 

3 Results. 

3.1 Lowest energy structure and isomers. 

In order to find the ground state structure we have to minimize the total en- 
ergy E of the cluster as a function of the positions Ri of all the ions. E is a 
function of 3iV — 6 independent variables (N is the number of ions) and min- 
imizing that function in a multidimensional space in which we must search 
for the absolute minimum among all the local minima is a tough problem. 
A restricted search can be performed by starting with a number of reason- 
able guesses for the geometry and allowing for the relaxation of those initial 
geometries. As input geometries we have considered cubic structures, rings 
and mixed structures, and for specific numbers of molecules we have also 
considered more open structures (like the wurtzite-type for n = 7, or a trun- 
cated octahedron for n = 12). All those structures are suggested by previous 
works using pair-potentials Jl2| or ab initio methods ||16||. The input geome- 



tries have been optimized with respect to a few parameters. More precisely, 
cuboid structures were relaxed with respect to a single parameter, the first- 
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neighbour cation-anion distance, whereas the ring structures were optimized 
with respect to two or three parameters: the cation and anion distances to 
the center of the ring, and the distance between rings. However, for clusters 
with 7 molecules or less, we have performed, in addition, a full relaxation of 
the structures obtained in the previous step. To this end, we have used a 
simplex algorithm p4j] , as given in reference [p^] . This is a simulated annealing 
recipe, which we have used with the "temperature" set to zero, resulting in 
the conventional downhill simplex method. The direct use of the simulated 
annealing technique would be prohibitive for our calculational resources. 

The results obtained for the geometries are shown in the figure 1. The 
Na + and I~ ions are represented by small and large spheres, respectively. 
The most stable structure (first isomer) for each size is shown on the left. 
The other structures are the low-lying isomers obtained (isomers several eV 
above the absolute minimum are not shown). In the case of clusters with 3, 
4 and 5 molecules we show two different views of the second isomer. The 
total energy difference with respect to the absolute minimum is given for 
each isomer. For n=2 we obtain a square as the ground state, for n=3 an 
hexagonal ring, for n=4 a cube, for n=5 a cube with a molecule attached 
to an edge and for n=6 an hexagonal prism. That is, we have found planar 
ground state structures for n = 2 and n = 3 and three-dimensional structures 
for n = 4 or larger. Planar ring-like isomers exist for n = 4 and n = 5, quite 
low in energy in the second case. The cuboid isomer for n = 6 is also quite 
close in energy to the hexagonal prism. The lowest energy structures for 
n = 2 — 6 are, however, not perfect. The full relaxation allowed by the 
calculation resulted in slightly distorted geometries. Nevertheless for n = 5 
the distortion of the cube is large: the cube edge facing the extra molecule 
enlarges its length. The structure for n = 5 can also be viewed as a distorted 
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(non planar) hexagon with two molecules attached, which can be justified in 
view of the structure of (Nal)^. The structure obtained for n = 7 is a piece 
of the wurtzite crystal. 

The same trend of competition between two different structural families 
continues for n larger than 7. n = 8, 11 and 14 are cuboid structures whereas 
n = 9, n = 10 and n — 15 correspond to stacks of hexagonal rings. Rather 
peculiar structures are obtained for n = 12 and n = 13: a hollow cage for 
n = 12 and a mixed structure (basket) for n = 13. The energy difference 
between the ground state and the cuboid isomers is so small for n = 10, 13 
and 15, that a more flexible relaxation of the cuboid isomer could perhaps 
change the order (we recall that the cuboid isomers have been relaxed with 
respect to a single parameter). This is unlikely for n = 9, for which the energy 
difference reported in Fig. 1 is 0.14 eV. In fact, we have performed additional 
calculations for the cuboid isomer in (Nal)g relaxing the energy with respect 
to more (six) parameters, but the hexagonal isomer still remained as the 
lowest energy structure. 



We proceed now to compare our results with pair potential calculations] 12 



In the pair potential calculations two different models have been considered: 
a rigid-ion model in which the energy is computed as a sum of the electro- 
static (Coulomb) and repulsive (Born-Mayer) energies, and a polarizable-ion 
model in which each ion has an induced dipole moment due to the fact that 
the closed electronic shells are polarized by the electric field of the others ions. 
We find discrepancies with respect to the rigid-ion model for n = 4, 5, 8, 11, 
14. In those cases the ground state predicted by that model corresponds to 
our second isomer. Those structures contain hexagonal or higher order rings. 
The same type of discrepancies are found compared to the polarizable-ion 
model for n — 8, 11, 14. Surprisingly, the polarizable-ion model favors the 
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defect-cuboid structure for n — 13, which is the second isomer in the PI cal- 
culations, only 0.02 eV above the ground state. This suggest the possibility 
of the defect-cuboid becoming the lowest energy structure in a PI calculation 
with a more flexible relaxation of interatomic distances. 

If we now compare the results of the PI model for (iVaC7) n ||16|| and (Nal) n 
we find a preference of the Nal clusters for more open structures, like the 
stacks of hexagonal rings for n — 6, 9, 10 and 15. Our calculations show 
the same general trends in ground state structures shown in other ab ini- 
tio calculations when the cation size is quite smaller than the anion size. 
For example, (LiF) n clusters studied in[^6j showed also a marked preference 
towards ring-like (open) structures. Thus, we can summarize these obser- 
vations in the following way: if the ions composing the cluster have similar 
sizes, or if the cations are larger than the anions, a preference for cubic struc- 
tures against hexagonal prismatic forms is advanced for those specific sizes in 
which both structures are feasible, as in the case of clusters with 6, 9 and 12 
molecules. On the other hand, if anions are larger than cations, a preference 
for hexagonal structures is advanced in the same cases. These hypothesis 
deserve further and more detailed studies on other ionic materials. 

The Nal crystal has the NaCl structure [|l|. The conclusion from the 
study of small (Nal) n clusters (n < 15) is that only for some values of n the 
cluster can be considered as a small piece of the bulk lattice. But there are 
other sizes, notably n — 6, 7, 9, 10, 12, 13, 15, for which the clusters, due to 
their very small sizes, prefer other structures. The magic numbers observed 
in the mass spectra of large Nal clusters []TU| (those clusters are, however, 
ionized and non-stoichiometric) have been explained as cuboids with the 
bulk structure. We conclude that the number of atoms in Nal clusters with 
15 molecules or less is not yet enough to produce bulk-like clusters. 
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3.2 Inter-ionic distances 

In figure 2 we present the average interionic distance d between nearest neigh- 
bours Na + and I~ ions for the ground state geometry as a function of the 
cluster size, d presents an irregular behaviour as a function of n, and this is 
due to the changes in the type of structure. By joining the points correspond- 
ing to the same type of structural family, either cubic or hexagonal, then the 
variation of d is less pronounced; this is shown by the two dashed lines. The 
cubic clusters have an interatomic distance larger than the hexagonal ones. 
d tends to a saturation value (~ 3.4A) that is higher than the distance in 
the bulk (db u ik(NaI) = 3.24A0). This can perhaps be attributed to the 
lack of electron correlation in the calculation. The incorporation of electron 
correlation tends to shorten the distances in ionic clusters, as other ab initio 



calculations 26] have shown for NaCl clusters. 



3.3 Binding energies 

The binding energy per molecule of the cluster with respect to the separate 
free ions, is given by: 

Eund = —[E(duster) - nE {r) - nE (Na + )] (13) 

where the energies of free Na + and Cl~ ions are - 161.67692 a.u. and - 
6918.06360 a.u, respectively. E bind is shown as a function of n in figure 3. 
The trend is an increase of the binding energy with n. But superposed to 
this general trend, especially stable clusters are predicted for n — (4), 6, 
9, 12. The stability of those clusters is reflected in local peaks in Eund (or 
a pronounced change of slope for n = 4). The magic numbers agree with 
those predicted earlier for (iVaCZ) ra using the PI model. From the num- 
bers given in figures 1 and 3 we notice that the energy differences between 
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isomers are very small compared to the total binding energies. This observa- 
tion is relevant because it indicates that the magic character of some clusters 
(the peaks in fig. 3) is associated to the specific number of molecules in those 
clusters and not to the particular structure of the lowest isomer. In fact, 
(NaCl) n clusters also have magic numbers for n = 6, 9, 12, although the 
lowest energy structures are cuboids (see fig.2 of ref|Tj|), that is, different 
from the lowest energy structures of (Nal) n . What makes special the num- 
ber of molecules in the magic clusters is that "compact" clusters can be built 
for those sizes. We illustrate this with specific examples. The groun state of 
(N 0,1)5 and the hexagonal isomers of (Nal)r contain some low coordinated 
ions, in contrast to (Nal)^. Both the ground state and the cuboid isomer of 
(Nal)g are more compact than the elongated forms of (Nal)g and (Nal) 10. 
Finally, the hexagonal and cubic isomers of (Nal) 12 are again more compact 
than the "defect" forms of (Nal)n and (Nal) 13. This idea of stability asso- 
ciated to compact clusters, which is evidently associated to the optimization 
of the attractive part of the electrostatic energy, provides justification for 
the interpretation of the periodic structure in the mass spectrum of large 
(NaI) n Na + clusters as revealing a cuboid pattern growth fL0| . 

3.4 Effects of the cluster size on the one-electron levels 
and vertical ionization potentials. 

In figures 4 and 5, the eigenvalues of the 2p orbitals of Na + and the 5p 
orbital of I~ ions are plotted as a function of n. The results correspond to 
the most stable structure for each cluster size. Instead of a single eigenvalue 
we have a band of eigenvalues for each cluster because the cations (or anions) 
occupy inequivalent positions in the cluster. In other words, the 2p (Na + ) 
and 5p (I~) eigenvalues are local atomic properties, and as such these could 
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determine the most reactive sites, or the preferred adsorption sites. 

The binding energy of the 2p level of the Na + ion, averaged over all the 
ions in the cluster, decreases slowly with cluster size and the average binding 
energy of the 5p level of the I~ increases. Of course, in the infinite crystal 
limit all the cations are equivalent, so there will be only one 2p eigenvalue 
for Na + (and a single eigenvalue for the 5p level of 

The clusters under study are made up of closed shells ions, whose wave 
functions are strongly localized. Thus, it can be assumed that an electron is 
extracted from a specific localized orbital when the cluster is ionized. This is 
the lowest bound 5p orbital, which generally corresponds to anions in corner 
positions. Using a Koopmans' like argument, we then identify the vertical 
ionization potential of the cluster with the smallest binding energy of a 5p 
electron. The dashed line in the figure 5 indicates the variation of the vertical 
ionization potential with the cluster size. 

The vertical ionization potential IP of a (Nal) n cluster in a pair potential 
model is a sum of two terms: the first is the energy necessary to remove an 
electron from the isolated anion, i.e., the electronic affinity EA of the neutral 
halogen atom, and the second is the electrostatic interaction of the electron 
removed, considered as a localized point charge, with the other ions in the 
cluster, also considered unit point-like charges, that is 

IP = EA + J2— (14) 

where is the distance between "i" and "j" sites. The electrostatic term 
provides a strong stabilizing effect (we note that the EA of a free I atom is 
only 3.06 eV^|). On the other hand, the short-range repulsive interaction 
among the electron and the electrons of neighbouring ions has been neglected 
in equation fll4|) . This interaction, which is taken into account in the PI 
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model, tends to lower the ionization energy below the value given by equation 

In order to study the different influence of the cluster field on cations 
and anions we show in Table I the binding energy differences A(nl) between 
the electronic levels (nl) of the ions in the cluster and in vacuum. As a 
representative cluster we have chosen the lowest isomer for n = 4, relaxed 
with respect the nearest-neighbour distance. The shifts A(ls), A(2s) and 
A(2p) are nearly the same for Na + , with differences in the meV range; this 
constant shift is the effect of the electrostatic term of equation (|i4|). However, 
the situation is different for the I~ anions, where the shifts are negative 
and larger for the outermost (5s and 5p) orbitals compared to the "core" 
orbitals. This is indicative of the high polarizability of the I~ anion. The 
distinct stabilizing effects experienced by the outermost orbitals of the I~ 
anions compared to "core" orbitals is ascribed to the quantum terms in the 
interaction energy, namely, the exchange interaction and the wave function 
orthogonality. 

3.5 Effects of the cluster formation on the wave func- 
tions. Partition of the binding energy 

The PI model assumes that the electronic density is a simple superposition 
of weakly overlapping spherically symmetric densities of the individual ions. 
Upon formation of the cluster, the Na + densities are practically equal to 
the free-ion densities, but a contraction of the I~ density is found which is 
mainly due to the 5p and 5s orbitals. This is in accordance with the general 
trends for ionic solids [17]. 

The orbital contraction of the outermost I~ orbitals increases with the 
cluster size, and is different for non equivalent anions. To get a more quan- 
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titative picture of this effect we proceed to study its influence on the cluster 
energy. The cluster binding energy of equation (13) can be written in a 
convenient form to our purposes using the equations (1) and (3): 

1 N 1 
E bind = --£ {Eti -E$+ -Ef nt ). (15) 

71 A=l Z 

By defining the deformation energy of the ion A as: 

Ei f = E* et - E&, (16) 
and the binding energy of the same ion as: 

Efind = -(Efef + 2 E f nt ), (17) 

we arrive to the following partition for the cluster binding energy: 

Ebind = -- Y,( E def + o E tt) = - E E f ind . (18) 
A Z n A 

According to this equation, we can separate the cluster binding energy into 
a sum of cationic and anionic contributions, each one composed in turn of 
two terms: one of them, E^, accounts for the energy associated to the de- 
formation of the ionic wave functions and the other is half of the interaction 
energy of that ion. The first term, E£,j, can never be negative: we just 
have to recall that E^ et and Eq are the expectation values of the intra-ionic 
part of the hamiltonian in states corresponding to the perturbed and the 
free ion, respectively. Since the last one minimizes that expectation value, 
we have E^ e t > Eq. Thus, the deformation energy is always an unfavourable 
contribution to the cohesion. In contrast, the interaction energy term, is neg- 
ative. The deformation and interaction parts of cationic and anionic binding 
energies are given in Table II for clusters with n = 1-9. P indicates the 
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specific type of site when there are several non equivalent positions. For 
n < 7, where we have made a "full optimization" of the structures, average 
values are quoted. We notice the passive cationic response to the cluster en- 
vironment, the corresponding deformation energies being negligible. Instead, 
there are sizable anionic deformation energies (up to 0.5 eV). The table dis- 
tinguishes differences between ions in vertex (V) and in edge (E) positions. 
The deformation energies are smaller for anions in vertex positions by 0.2 
eV. This is related to the lower coordination for ions in vertex positions. 
However, in all cases, the anionic interaction energies are larger, in absolute 
value, than the corresponding deformation energies and their sum is always 
negative. Another interesting aspect is that we have at our disposal the 
binding energies of the ions at all nonequivalent positions in the cluster. 

The binding energies in figure 3 grow quite fast with n, but the limited 
size range studied here is not enough to estimate precisely the limiting value. 
That limit seems to be a little smaller than the experimental binding en- 
ergy in the bulk crystal (7.05 eV/molecule[l|]). Our theoretical value for the 
molecule is also smaller than the experimental one. Three factors can con- 
tribute to those errors. In the first place, the geometries for n > 7 have 
not been fully optimized. Second, there is an intrinsic variational restric- 
tion in the PI method, which is associated to the neglect of covalent mixing 
between anionic and cationic wave functions. And third, we have neglected 
the coulomb electronic correlation. We think that the main source of error 
is the last one, in line with the arguments leading to the same conclusion in 
alkali halide solids|T7j]. Although it is sometimes stated that the HF approx- 
imation would give good binding energies for ionic compounds because the 
correlation errors in the free-ion and cluster wave functions tend to cancel 
out, the orbital anionic deformation upon cluster formation casts doubts on 
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this assumption. 

4 Conclusions 

(Nal) n clusters up to n = 15 have been studied using the ab initio Perturbed 
Ion (PI) model. Geometries have been "fully" relaxed only for n < 7, due to 
the increasing computational time necessary to carry out these optimizations 
in larger clusters. The isomeric study shows that these clusters are not 
yet fragments of the bulk crystal lattice, although lattice fragments occur 
indeed for some values of n. Comparison to pair potential models shows some 
discrepancies. Interionic distances do not show a smooth variation unless 
we restrict ourselves to a fixed geometrical family. Though extrapolation 
to very large cluster sizes is difficult, we think that inclusion of correlation 
effects would be important in order to obtain quantitative agreement with 
bulk quantities. Particularly stable ("magic") clusters have been obtained 
for n = (4), 6, 9 and 12. 

Cationic electron densities remain nearly undisturbed upon formation of 
the cluster. On the contrary, anionic densities are contracted in this process. 
The orbitals affected have been identified as the outermost ones, mainly 
the 5p orbital of I~ . Nevertheless the 5s orbital also experiments some 
deformation, and should not be considered as a "core" orbital in a self- 
consistent-field calculation. This effect has been analyzed by performing a 
study of the energy eigenvalues. In the PI model, the binding mechanism 
is a combination of the orbital deformation and the quantum mechanical 
interaction between each ion and the cluster, and the cluster binding energies 
have been analyzed from the point of view of this bonding mechanism. 
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Captions of tables 



Table I. Differences among the binding energies of the orbitals in (N 0,1)4 
and the respective values for the ions in vacuum. The shifts are given in eV. 

Table II. Cation and anion deformation, interaction, and binding ener- 
gies for several cluster sizes (in eV). The position P is given explicitly if 
there are different types of sites: edge (E), vertex (V). 
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Captions of figures 



Figure 1. Isomeric geometries for (Nal) n clusters. The most stable 
structure (first isomer) is shown on the left side. Total energy differences 
(in eV) with respect to the most stable structure are given for each isomer. 
Na + , small spheres; I~, large spheres. For n = 3 — 5, two different views of 
the second isomer are provided. 

Figure 2. Averaged nearest neighbour Na-I distance for the structures 
of lowest total energy in Figure 1. The two lines join cubic and hexagonal 
clusters, respectively, n = 7, 12, 13 correspond to peculiar structures. 

Figure 3. Binding energy per molecule as a function of the cluster size. 

Figure 4. Orbital energies (with opposite sign) of the 2p levels of Na + 
cations as a function of the cluster size for the ground structures of Figure 1. 

Figure 5. Orbital energies (with opposite sign) of the 5p levels of I an- 
ions as a function of cluster size for the ground structures of (Nal) n clusters. 
The dashed line joins the vertical cluster ionization potentials. 
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Na+ 


(ls) 
+6.239 


(2s) 
+6.237 

(2 P ) 
+6.232 








r 


(Is) 


(2s) 


(3s) 


(4s) 


(5s) 




-4.697 


-4.714 


-4.728 


-4.758 


-5.251 






(2p) 


(3p) 


(4p) 


(5p) 






-4.714 


-4.731 


-4.766 


-5.241 








(3d) 


(4d) 










-4.725 


-4.768 





Table 1: Table I. 
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cations 


anions 


n 






^bind 


P 


Erfe/ 




^bind 


P 


1 


0.0003 


- 4.8238 


2.4116 




0.1301 


- 3.8487 


1.7942 




2 


0.0002 


- 5.6798 


2.8397 




0.2211 


- 5.0630 


2.3104 




3 


0.0001 


- 6.0880 


3.0439 




0.2582 


- 5.2000 


2.3418 




4 


0.0002 


- 6.1880 


3.0938 




0.3612 


- 5.7750 


2.5263 




5 


0.0002 


- 6.3740 


3.1868 




0.3217 


- 5.7310 


2.5438 




6 


0.0005 


- 6.4750 


3.2370 




0.3620 


- 5.8300 


2.5530 




7 


0.0008 


- 6.4180 


3.2082 




0.3684 


- 5.7100 


2.4866 




8 


0.0004 


- 6.3380 


3.1686 


V 


0.2935 


- 5.7570 


2.5850 


V 




0.0003 


- 6.5350 


3.2672 


E 


0.4861 


- 6.0110 


2.5194 


E 


9 


0.0001 


- 5.9640 


3.1915 


V 


0.3863 


- 4.6450 


2.4652 


V 




0.0001 


- 6.5060 


3.4721 


E 


0.5969 


- 7.8230 


2.5505 


E 



Table 2: Table II. 
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